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Abstract 

We study the effects of tlie tidal interaction with the companion, via orbital separation and binary 
mass ratio, on the global one-armed oscillation modes in disks around binary Be stars. Our model takes 
into account the three-dimensional effect that contributes to the mode confinement, which was recently 
found by Ogilvie (2008). We find that the one-armed oscillations are well confined in systems with disks 
larger than a few tens of stellar radii. In such systems, the oscillation period depends little on the binary 
parameters. On the other hand, in systems with smaller disks, where the mode confinement is incomplete, 
the oscillation period increases with increasing orbital separation and/or decreasing binary mass ratio. The 
eigenmode is insensitive to the spectral type of the central star. Our results suggest that the dependence 
of V/R oscillation period on the orbital separation and binary mass ratio should be observed only in short 
period binary systems, and that, for systems with a similar orbital period, those with higher mass ratios 
will show shorter V/R variations. 
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1. Introduction 

Be stars are non-supergiant B-type stars whose spectra 
show, or have at some time shown, one or more Balmer 
lines in emission (Collins 1987). A Be star is a rapidly 
rotating, B-type star as the central star with a two- 
component circumstellar envelope, a polar wind and an 
equatorial disk. The polar wind is a low-density, fast out- 
flow emitting UV radiation. The wind structure is well ex- 
plained by the line-driven wind model (Castor et al. 1975; 
Friend & Abbott 1986). On the other hand, the equato- 
rial disk is a geometrically thin, high-density plasma in 
nearly Keplcrian rotation (e.g.. Porter & Rivinius 2003), 
from which the optical emission lines and the IR excess 
arise. Although there are still several competing scenarios 
for the Be disk formation, the most promising one is the 
scenario where the disk is formed by viscous decretion of 
gas ejected from the central star (Lee et al. 1991; see also 
Porter & Rivinius 2003, and references therein). 

Many Be stars show long-term variations in their spec- 
tra over years to decades. One of these phenomena is 
called long-term V/R variations, which are variations of 
the ratio of relative intensity of violet (V) and red (R) 
peaks of a double-peaked emission line profile. The pe- 
riod of the long-term V/R variations is typically in the 
range 5-10 yr for isolated Be stars. It is widely accepted 
that the long-term V/R variations are attributed to global 
one-armed (i.e., the azimuthal wave number m — 1) oscil- 
lations in the Be disk (e.g.. Porter & Rivinius 2003). It 
is based on the fact that in a nearly Keplerian disk a 
one-armed mode shows up as a very slowly revolving per- 



turbation pattern (e.g., Kato 1983). The early versions 
of the one-armed oscillation model (Okazaki 1991, 1997; 
Papaloizou et al. 1992; Savonije & Heemskerk 1993) qual- 
itatively well explained the observed characteristics of the 
V/R variations. However, owing to the lack of the mech- 
anism for confining the modes to the inner part of the 
disk and the high sensitivity of the oscillation period on 
various parameters, the model had little predictive power 
(Fift & Harmanec 2006). 

Recently, the model was advanced greatly. Ogilvie 
(2008) has solved the mode confinement problem by in- 
vestigating a three-dimensional effect on the mode char- 
acteristics. While previous versions of the model had as- 
sumed motions to be independent of height since the disk 
is geometrically thin, he argued that the variation of the 
vertical gravitational acceleration around an elliptical or- 
bit excites an oscillatory vertical motion in an eccentric 
disk that should not be neglected. Ogilvie (2008) showed 
that the three-dimensional effect alone allows confined 
prograde modes. 

Hitherto, theoretical studies of one-armed oscillation 
modes in Be disks are mostly limited to isolated Be stars. 
However, according to the current census, the fraction 
of Be stars found in binary systems is ~ 1/3 (Porter & 
Rivinius 2003). The presence of the companion is most 
likely to change some features of the global one-armed 
oscillations in Be disks. Actually, in Be/X-ray binaries, 
which consist of a neutron star and an early-type Be star, 
the time-scale of V/R variations is frequently about lyr, 
which is much shorter than the typical V/R time-scale for 
isolated Be stars (e.g., Reig et al. 2005). 
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Moreover, there are observationally two types of long- 
term V/R variations in binary Be stars (Stefl et al. 2007). 
One is quasi-periodic V/R variations, where V/R peak 
ratio varies quasi periodically, which shows no correlation 
with their orbital phase. This type of V/R variations are 
similar to those usually observed for V/R variations in iso- 
lated Be stars. The other type is periodic V/R variations 
locked to the orbital period. This characteristic is found 
only in binary Be star systems. Unfortunately, no the- 
oretical explanation is available yet for the phase-locked 
V/R variations. 

In this paper, we investigate the effects of the compan- 
ion on the characteristics of one-armed oscillation modes 
in binary Be stars, taking into account the tidal effects in 
addition to the three-dimensional effect. For simplicity, 
we assume the binary orbit to be circular. We find that 
the mode is well confined in disks larger than a few tens 
of stellar radii, which is consistent with Ogilvie (2008). In 
smaller disks, however, the mode confinement is incom- 
plete and the oscillation period depends on the binary 
parameters significantly. 

The structure of this paper is as follows. In section 
2, we describe the unperturbed disk and derive the basic 
equations for linear, one-armed perturbations, following 
the formulation by Ogilvie (2008). The boundary condi- 
tions to be adopted are also discussed. In section 3, we 
present our numerical results, and in section 4, we com- 
pare our results with the observational features of V/R 
variations in binary Be stars. The final section is devoted 
to conclusions. 

2. Basic equations for m — 1 density waves 

2.1. Basic equations 

As an unperturbed state, we take a geometrically 
thin, axisymmetric disk that locally rotates at a nearly 
Keplerian speed and is in hydrostatic equilibrium in the 
vertical direction. Because of the strong photoionization 
heating from the stellar radiation, this Be disk is taken 
to be isothermal at O.BTeff (Carciofi & Bjorkman 2006), 
where T^s is the effective temperature of the central star. 
For simplicity, we assume the binary orbit to be circular 
and neglect the advective motion and viscous effects in 
the Be disk. 

We use the cylindrical coordinate system {r,<f>,z). The 
equation of continuity and the equation of motion are re- 
spectively given by 



| + v.(H^o 



and 



(1) 



(2) 



Here v is the velocity vector, ^E" is the gravitational po- 
tential, and h = In p is the enthalpy of an isothermal 
gas, where Cg is the isothermal sound speed and p is the 
density. 

As far as eigenmode oscillations are concerned, we may 



consider only the azimuthally averaged tidal potential 
(Hirose & Osaki 1993). We also take into account the 
rotational deformation of the rapidly rotating Be star, in- 
cluding only monopolc and quadrupole term (Papaloizou 
et al. 1992). The potential in the disk midplane, 'I'„j, is 
then reduced to 
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where Mi, Ri and f2i are the mass, radius an d angular 
rotation speed of the Be star, k2 and flc = a/ GMi/R\ 
are its apsidal motion constant and critical angular rota- 
tion speed, M2 is the mass of the companion, and D is 
the binary separation. In equation (3), the first term is 
the gravitational potential of the Be star, in which the 
quadrupole contribution due to the rotational distortion 
of the star is taken into account, and the second term is 
the azimuthally averaged tidal potential. 



In circular binaries with small mass ratios {M2/M1 



< 



0.3), a viscous disk is truncated at the tidal radius given 
by 



R, 



tides ' 
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(4) 



(Whitehurst & King 1991). Here Rl is the Roche lobe 
radius of the Be star given approximately by 
0.49g-2/3 
^'^"■^0.69g-2/3 + ln(l + g-i/3) 

(Eggleton 1983); where q ~ M2/M1 is the binary mass 
ratio. In this paper, we thus assume that the Be disk is 
truncated at r = i?tidcs- 

The unperturbed equilibrium state of the isothermal 
disk is given by 

-v2 



P = Pm(f)exp 



2m 



t;=[0,rf](r),0]. 



and 
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(6) 

(7) 
(8) 

(9) 



dr dr 

where fi(r) is the angular frequency of disk rotation, Pm{f) 
is the midplane density, and H{r) = Cs/i^K with f^K = 
\J GMi I is the scale-height of the disk. 

From equations (3) and (9), we have the explicit form 
of as follows. 

1/2 



(10) 



where / = ili/ilc. Then, the associated local epicyclic 
frequency K(r) is explicitly written as 
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2n (2n + r^ 
dr 

m 
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(11) 



On the above unperturbed state, we superpose a lin- 
ear, m = 1 isothermal perturbation in the form of normal 
mode of frequency w, which varies as exp[z((/) — ^jJt)]. The 
linearized perturbed equations arc then obtained as fol- 
lows. 
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(12) 
(13) 
(14) 

(15) 



where {v'^^v'^^v'^) and h' are the perturbed velocity and 
enthalpy, respectively. 

In order to take the three-dimensional effect into ac- 
count, we expand perturbed quantities in the z-direction 
in terms of Hcrmitc polynomials as 



v'r{T,z) = ^U„(r)if„(C), 

n 

v'^{r,z)=Y,yn{r)H,,{C), 

71 

n 

/i'(r,z)=^/i„(r)if„(C) 



(16) 
(17) 
(18) 
(19) 



(Ogilvie 2008; see also Okazaki et al. 1987), where Hn{C,) 
is the Hermite polynomial defined by 



HniO 



exply 



(20) 



with ^ = z/H being a dimcnsionlcss vertical coordinate 
and n = 0, 1,2, . . .. 

The system of resulting equations is not closed because 
the equation for u„ refers to the higher order term hn+2, 
which depends on Un+2, the equation for which in turn 
refers to /i„+4, and so forth. Following Ogilvie (2008), we 
close the system of resulting equations by assuming u„ = 
for ri > 2. Neglecting U2 compared to Uq is equaivalent to 
assuming that the eccentricity of perturbed orbit of each 
gas particle is independent of z. If /12 is also neglected, 
then we will obtain the two-dimensional equations. The 



basic equations for linear m ~ 1 perturbations in inviscid 
disks arc then given as follows. 
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(21) 
(22) 
(23) 



where E = (27r)^/^p„iiJ is the surface density and is 
the local apsidal precession frequency given by 




(24) 



In deriving equations (21)-(23), we have used approxima- 
tions |w| < Jl, |ri - k| < fl, and (cs/rO)^ < 1. Note that 
these equations are essentially the same as equations (A6)- 
(A8) of Ogilvie (2008), except that our equations implic- 
itly depend on the tidal effect through and k. 

Eliminating h2 from equations (21) and (23), we have 

9c? "I 

2 (w — cijp ' 
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dr 
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and equation (22) is written as 
.duQ Q (1 dhiT,\ luo 



rfln?' 



We define Yi and Y2 as Yi = Mq and Y2 = - 
the basic equations to be solved are given by 



dr 

and 
dY2 
dr 



f3 



Yi 



,Y2 



1±\ 

Ar'^Vl ) 



Yi - -I2, 
r 



(25) 

(26) 
-iho- Then 

(27) 
(28) 



where /3 = — 1/2 — (dlnS/dlnr). Note that as shown by 
Ogilvie (2008), the term 9Cg/4r^ri results from the three- 
dimensional effect and provides an important contribution 
to the confinement of the m = 1 modes. 

2.2. Boundary conditions 

We now consider the boundary conditions for equations 
(27) and (28). In Be stars with growing or persistent 
equatorial disks, material is likely injected continuously 
from the star to the disk, so there is no gap between the 
star and the disk. Since the pressure scale-height of the 
Be star near the interface is much smaller than that of 
the disk, the oscillations cannot penetrate the stellar sur- 
face. Hence, throughout this paper, we take the rigid wall 
boundary condition at the disk inner radius, i.e., u'. = 
at r ~ Ri, except in section 3.4 where we study the effect 
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of the free inner boundary condition. With Yi — uq, this 
condition is written as 

Yi^Q at r = Ri. (29) 

In our model the Be disk is truncated at the tidal radius. 
Hence, we take the free boundary condition at the outer 
boundary. This condition is written as 

Ap = 0, (30) 

where Ap is the Lagrangian perturbation of pressure. 
Given that the Eulerian perturbation of pressure, p' , 
varies as exp[z(0 — wt)], equation (30) can be written as 

^(^7-c.)y + ^;;^ + ^;;^=0. (31) 
or oz 

Using the same procedure as in the previous subsection, 
we have 

iin-iu)ho + uo-^ + uoc^^— = 0, (32) 
p or or 

where we have used h' — p/ p. Finally, with Yi — uq and 
Y2 = —iho, we have the outer boundary condition as 

dlnE 

{n~uj)Y2-- ^^1=0 at r = i?tides. (33) 

dlnr r 

3. Numerical results 

We solve basic equations (27) and (28) for a range of 
parameters, including binary separation D (10-Ri < D < 
50i?i) and binary mass ratio q {q — 0.1 and 0.3), with 
the boundary conditions described above. In the follow- 
ing calculations, we take a BOV star with Mi = IT.SA/q, 
Ri = 7.4i?0, and Tee = 30000K (Cox 2000) as the Be 
star model, except in section 3.3 where we consider a 
B5V star to study the spectral dependence of the mode 
characteristics. As for the quadrupole parameter, we 
take ^2/^ = 5 X 10"'^ as a representative value, although 
there is a large uncertainty of this parameter in the range 
^2/^ = 2 X 10"'^ — 10~^, as discussed in the next section. 

The pressure gradient force in our isothermal disk model 
depends solely on the radial density distribution. In the 
following we take the power-law disk density distribution 
that varies as pm oc r~^/^, which is the theoritical density 
distribution of steady, isothermal decretion disks (Porter 
1999; Okazaki 2001; Carciofi & Bjorkman 2006). 

In this section, we study the effect of orbital separation 
in section 3.1 and that of binary mass ratio in section 3.2. 
We then briefly discuss the results for B5V stars in section 
3.3 and the effect of the free inner boundary condition in 
section 3.4. 

3.1. Effect of binary separation 

In order to study the effect of binary separation on the 
oscillation characteristics of the Be disk in circular bina- 
ries, we vary the value of the binary separation D from 
D = 10i?i to D = 50i?i, with binary mass ratio q = 0.1. 

Figures 1 and 2 show the fundamental mode (left) and 
the first overtone (right) for D — 10i?i and D = 50i?i, re- 
spectively. The top panels in each flgure show the cigen- 
function Yi (= uq) (solid hne) and Y2/c^ (= —iho/cl = 



[Vol. , 

—ipo/p) (dashed line) in the radial direction. In the mid- 
dle panels, the color-scale plot shows the surface density 
perturbation relative to the unperturbed surface density, 
E'/S, in the (r, 0)-plane: the red (blue) region has a pos- 
itive (negative) density perturbation. The arrows super- 
posed on the color-scale plot denote the velocity vectors 
associated with the mode. The disk itself rotates counter- 
clockwise. The bottom panels shows the density perturba- 
tion, p' , in the (r, z)-plane of the disk, at (/> = 0. Note that 
the modes are linear, so their normalization is arbitrary. 

All modes in Figs. 1 and 2 have positive frequencies, 
that is, they are prograde modes, which precess in the di- 
rection of disk rotation. Figure 2 (left) shows that the fun- 
damental mode is well confined in a large disk. However, 
as shown in Fig. 1 for D = 10i?i, when the binary sepa- 
ration is small, the confinement of the fundamental mode 
is weak, so the mode propagates over the whole disk. As 
for the first overtones, the confinement is weak: it is not 
confined well even for a wide binary with D ~ bORi . Note 
that these results are consistent with Ogilvie (2008) who 
showed that the 3-D theory allows prograde, fundamental 
modes to be confined within several tens of stellar radii in 
disks around isolated Be stars. 

In the middle panels of Figs. 1 and 2, we note that per- 
turbation pattern averaged in the vertical direction is sim- 
ilar to that for 2-D disk models around isolated Be stars, 
which have been extensively studied. We thus expect that 
the current model can also explain the observed line profile 
variability. In fact, the modal features suggest that when 
the violet (red) peak is stronger than the red (violet) peak, 
the line profile as a whole redshifts (blueshifts). This is 
a typical behavior of the observed line profile variability 
(e.g., Okazaki 1991). 

The bottom panels of Figs. 1 and 2 show a new feature 
that has not been suggested by the previous 2-D stud- 
ies: The density perturbation, p'{r,z), does not always 
take the maximum in the disk midplane. The behavior of 
p'{r,z) depends on that of the radial component of per- 
turbed velocity, wo(r), as seen from eqauation (23). In 
regions where ug is negligible, p' is dominated by the con- 
tribution from Hq, so\ p' \ 'At a fixed radius has a maximum 
in the equatorial plane and decreases with vertical coor- 
dinate. In contrast, in regions where ito is not negligible, 
the contribution from h2 becomes significant. Then, the 
density perturbation has a maximum at a vertical coordi- 
nate z ~ (1 — 2)H , not in the equatorial plane. In the first 
overtone, these patterns apear alternately in the radial 
direction. 

The distribution of eigenfrequency is shown in Figure 3 
as a function of binary separation normalized by the Be 
star radius, D/Ri, for q = 0.1. Here and hereafter, we 
normalize the eigenfrequency to by the stellar critical ro- 
tation frequency fie = \/ GMi/ R\. The solid lines and 
the dashed lines denote the fundamental modes and the 
first overtones, respectively. The thick lines are for a BOV 
star (and the thin lines are for a B5V star discussed in 
section 3.3). Note that all modes arc prograde modes. 
Note also that the oscillation period of the fundamental 
mode increases to an asymptotic value with increasing bi- 



F. Oktariani and A. T. Okazaki 



No. 



Global Disk Oscillations in Binary Be Stars 



>- 







\ 





\ 

\ 




: 


1 f 1 1 1 1 1 1 1 1 

\ 

\ 




: \ 




\ 




. Fundamental 

.... 1 .... 1 


■node 





-4 -2 


2 4 




-4 -2 


2 4 


x/R, 






x/R, 














1 -0.5 


0.5 

r/z 


- 


1 -0.5 


0.5 






3 






1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 


1 1 1 1 1 1 1 1 1 


'''I'll 






2 

X 

\ 

N 

1 






i 


' 









2 3 


4 


2 3 


4 


r/R, 






r/R, 














1 -0.5 


0.5 

P 


- 


1 -0.5 


0.5 

P 



Fig. 1. Eigcnfunctions of the m = 1 fundamental mode (top 
left) and the first overtone (top right) for D = lOiJi and 
q = 0.1 with the rigid inner boundary condition. The solid 
line denotes Yi, while the dashed line denotes ¥2/0^. In 
the middle panels, the color-scale plot show the surface 
density perturbation normalized by the unperturbed surface 
density, T,' /T,, in the (r, i;i)-plane, while the arrows denote 
the perturbed velocity vectors. The bottom panels show 
the density perturbation, p', in the (r, 2)-planc at = 0. 
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Fig. 2. Same as Fig. 1, but for D = 50Ri. 

nary separation: for a BOV star, it varies from 1.38 yr for 
D = 10i?i to 2.25 yr for D 50i?i. This is because the tidal 
field is weaker for larger orbital separations, and weaker 
tidal field decreases the local apsidal precession rate uup^ 
[see eq. (24)]. The lower the local precession rate uup^, the 
lower the eigenfrequcncy uj. Thus, the cigcnfrequcncy of 
the fundamental mode decreases with increasing orbital 
separation. However, the effect of the binary separation 
is only appreciable for small D /Ri . In systems with large 
binary separation, the eigenfrequcncy of the fundamental 
mode changes little, because the mode is already well con- 
fined to the inner part of the disk, so the size of the disk 
does not matter anymore. 

3.2. Effect of binary mass ratio 

In order to study the effect of binary mass ratio on disk 
oscillation modes, we have performed calculations for q = 
0.3, in addition to the calculations for q = 0.1 shown above. 
We have also carried out calculations for an artificial case 
with no tidal effect, except that the disk is truncated at 
the same radius as for q = Q.l (hereafter, called the q = 
case), to compare with the q = 0.1 result. 

Figure 4 shows the distribution of eigenfrequcncy as a 
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Fig. 3. Eigenfrequency as a function of normalized binary 
separation, D/Ri, for q = 0.1. The solid linos denote the 
fundamental modes, while the dashed lines denote the first 
overtones. The thiek lines and the thin lines are for a BOV 
star and a B5V star, respectively, ui is normalized by the 
stellar critical rotation frequency Qc. The rigid boundary 
condition is applied at the inner disk radius. 
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Fig. 4. Eigenfrequency as a function of normalized binary 
separation, D/Ri, for q = (thin lines), 0.1 (lines with in- 
termediate thickness), and 0.3 (thick lines). The solid lines 
are for the fundamental modes, while the dashed lines for the 
first overtones, lj is normalized by the stellar critical rotation 
frequency Qc. The rigid inner boundary condition is adopted. 



function of normalized binary separation, D / Ri , for q~0 
(thin lines), 0.1 (lines with intermediate thickness), and 
0.3 (thick lines). The solid line and the dashed line de- 
note the fundamental mode and the first overtone, respec- 
tively. It is noted that the eigenfrequency increases with 
increasing binary mass ratio. This is particularly the case 
for small D/Ri. This is because the local apsidal preces- 
sion rate Wpr increases with increasing mass ratio q. As 
mentioned above, the higher Wpr, the higher the eigenfre- 
quency. 

Here we compare our result with that of Ogilvie 
(2008). Ogilvie (2008) normalized the eigenfrequency lo 
and the quadrupole parameter ^2/^ by , where e = 
Cs(i?i/GMi)V2^ which is ~ 2.3 x lO'^ in our case of the 
BOV star. Then, fca/^ = 5 x lO^^ and w/f2c 6.8 x 10"'' 
for q = {) and D = 50i?i in our case corresponds to 
Q = k2p/e'^ ~ 9.1 and uj = Lu/e'^ ~ 1.3, respectively. We 
note that this is consistent with a) ~ 1.2 for Q ~ 9 as seen 
by eye from Fig. 3 (left) of Ogilvie (2008)). 

Figure 4 suggests that the first overtones are retro- 
grade modes, which precess in the opposite direction to 
disk rotation, for small binary separation {D ~ 10i?i) and 
very small mass ratio {q <Si 0.1). Note that the retrograde 
modes can exist in the Be disk in binaries because of its 
truncated outer radius. 

3.3. Spectral type dependence 

As shown by Fift & Harmanec (2006), the one-armed 
oscillation period depend also on the mass and radius of 
the central star. In this subsection, we discuss the effect 
of the spectral type. For this purpose, we calculated the 
eigenmdoes for a B5V central star with Mi = 5.9Mq, Ri = 
3.9i?o, and T^s = 15200K (Cox 2000). 

Figure 3 compares the distribution of normalized eigen- 
frequency for two different spectral types, BOV (thick 
lines) and B5V (thin lines). From the figure, we note that 
the normalized eigenfrequency of the fundamental mode 
for a B5V star is slightly lower than that of a BOV star. 
This is because both the local apsidal precession rate LUp^ 
and the 3-D effect are smaller for a B5V star than for a 
BOV star: A B5V star has a higher critical rotation fre- 
quency fie and a lower sound speed than a BOV star does. 
This decreases ojpr/flc- It also weakens the 3-D effect via 
a smaller value of 9Cg/4r^J7 in eq. (28). The same dis- 
cussion is applied to the first overtone, except for small 
binary separation {D ^ 20i?i) where the tidal effect for a 
fixed value of q is relatively stronger for a B5V star than 
for a BOV star. 

Because of a higher ftc of a B5V star, the dimensional 
eigenfrequency (or the oscillation period) shows the op- 
posite trend. Although the normalized eigenfrequency is 
slightly lower for a B5V star than for a BOV star, the 
eigenfrequency itself shows the opposite trend. Oscillation 
period of the fundamental mode for a B5V star ranges 
from 0.96 yr for D = 10i?i to 1.5yr for D = 50i?i, which 
is slightly shorter than that for a BOV star by ~ 30%. 
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Fig. 5. Eigenfrequency as a function of normalized binary 
separation, D/Ri, for free (thick lines) and rigid (thin lines) 
inner boundary condition. The solid lines are for the funda- 
mental modes, while the dashed lines for the first overtones. 

3.4- Free inner boundary condition 

In long term, some Be stars change their state between 
the Be star phase and the normal B star phase. There 
are also other stars that show the loss and reformation of 
the high-velocity wings of emission lines (Rivinius et al. 
2001). These observational features indicate that in these 
stars, the Be disk, or at least the inner part, is completely 
lost and then reformed (Porter & Rivinius 2003). In the 
formation stage, it is likely that the Be disk and the cen- 
tral star are directly connected. The rigid inner boundary 
condition is adequate for such a situation. In the dissi- 
pation stage, however, a gap is expected to open between 
the disk and the central star, owing to the accretion or 
ablation of the innermost part. 

In order to study global oscillations in dissipating Be 
disks, we have performed calculations for the free inner 
boundary condition, which means that the Lagrangian 
perturbation of the pressure, Ap, vanishes at the disk in- 
ner radius. For simplicity, we have adopted the same un- 
perturbed disk state as in the rigid boundary case. Using 
the same procedure as in section 2.2, the free inner bound- 
ary condition is written as 

{^-uj)Y^-— ^yi==0 at r = i?i. (34) 

rfmr r 

Figure 5 shows the eigenfrequency distributions as a 
function of normalized binary separation, D / Ri^ for q = 
0.1. The thick lines are for the free inner boundary con- 
dition, while the thin lines for the rigid inner boundary 
condition. The solid lines and the dashed lines denote the 
fundamental modes and the first overtones, respectively. 
We note that the eigenfrequency is higher and the mode 
is more concentrated in the inner part with the free in- 
ner boundary condition than with the rigid one, as first 
pointed out by Papaloizou & Savonije (2006) and then 
confirmed by Ogilvie (2008). As in the case of rigid inner 



boundary condition, our result for q = 0.0 and D — 50i?i 
with the free inner boundary condition is consistent with 
that of Ogilvie (2008). 

4. Discussion 

As mentioned in section 1 , binary Be star systems show 
two types of long-term V/R variations. The first one is 
quasi-periodic V/R variations, similar to those observed 
in isolated Be stars. The second one is periodic V/R vari- 
ations locked to their orbital period (Stefl et al. 2007). 
In our study, we have found only eigenmodes compatible 
with the former type of variations: our eigenmodes de- 
pend on the orbital period, but they are not locked to 
it. It is interesting to carry out numerical simulations to 
study whether tidally forced, nonlinear oscillations can ex- 
plain the phase locked V/R variations, but this is beyond 
the scope of this paper. 

In the previous section, we have shown how global one- 
armed oscillation modes depend on the binary parame- 
ters. The oscillation period becomes longer in the case of 
a wider binary separation and/or a lower binary mass ra- 
tio. Such dependence, however, will be observed only for 
short period binaries. In wide binaries where the disk size 
is large enough to confine the one-armed modes in the 
inner part of the disk, the mode characteristics depend 
little on the binary parameters. The critical orbital sepa- 
ration, where the size of the Be disk is just large enough 
to confine the fundamental mode, depends on the details 
of the model. The difference between the model oscilla- 
tion periods (a few yr) and the observed ones (5-10 yr) 
implies that there are still missing mechanisms that affect 
the global disk oscillations. One of them could be the op- 
tically thick line forces examined by Gayley et al. (2001), 
who found that these forces have a net effect that likely 
lowers the eigenfrequency of the one-armed mode. The 
more realistic density and temperature distributions could 
also contribute to lower the eigenfrequency. However, de- 
spite these mechanisms not taken into account, we believe 
that the current results are qualitatively robust. 

Observationally, no clear dependence of the V/R pe- 
riod on the binary parameters has been found. Reig 
et al. (2005) compared the time-scales of disk forma- 
tion/dissipation cycle and V/R variability with the or- 
bital period for eight Be/X-ray binaries (see their Table 3). 
They found a good correlation between the time-scale of 
disk formation/dissipation cycle and the orbital period, 
which agrees with the scenario that the Be disk is trun- 
cated by the tidal torques from the neutron star (Okazaki 
& Negueruela 2001). On the other hand, the V/R vari- 
ability time-scale has no clear correlation with the orbital 
period, although short period systems have, in general, 
short V/R time-scale, with which our results qualitatively 
agree. 

In this paper, we have assumed an inviscid disk for 
mathematical simplicity. However, viscosity in the Be disk 
is considered to play an important role in the formation 
of a nearly Kcplerian disk by causing the outward drift of 
matter ejected from the equatorial surface of the Be star. 
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as proposed in the viscous decretion disk model by Lee et 
al. (1991). Negueruela et al. (2001) showed, after Kato et 
al. (1978), that global m = 1 modes, which arc neutral in 
inviscid disks, become ovcrstablc when the viscous effect 
is taken into account as a perturbation. The growth rate 
of the mode, i.e., the imaginary part of the eigenfrequency, 
of the mode is proportional to the Shakura-Sunyaev vis- 
cosity parameter a, while the real part is independent of 
a as long as a ^ 1. Ogilvie (2008) also obtained the result 
that the real part of an eigenfrequency in a viscous disk 
with a = 0.1 is only slightly different from that in an in- 
viscid disk. Thus, we expect that as far as the oscillation 
period is concerned, our result is valid for viscous disks 
with small a as well. 

Finally, we briefly comment on the effect of stellar ro- 
tation via the quadrupole parameter ^2/^. For B-type 
main sequence stars, the stellar evolution calculations 
result in the apsidal motion constant ^2 in the range 
2.5 X 10"'^ < ^2 < 10~^, depending on the evolutionary 
stage and the internal angular momentum distribution 
(Claret & Gimnez 1991; Claret 1995; Claret 1999). While 
it remains controversial how close the rotation of Be stars 
is to the critical rotation, it is likely / = Oi/57c ^ 0.9 (e.g., 
Fremat et al. 2005), which yields for the quadrupole factor 
in the range 2 x lO'^ < < 10"^. 

In order to study the effect, we have calculated the 
eigenmodes for ^2/^ = 10^^, implicitly assuming an ex- 
tremely high value of fc2 ^ 10^^ for a critically rotating 
star (/ ~ 1), for a BOV star in a binary with q = 0.1. As 
expected, the resulting eigenfrequency was higher than 
that for = 5 x 10~^: It increased from 1.1 x 10~^ 
to w ~ 1.6 X 10-^ for D = 10i?i and from w ~ 6.8 x lO"'^ 
to 1.3 X 10"'^ for D ~ 50i?i. This is because the local 
apsidal precession rate LUp,- apprcciablly increases in the 
inner part with an increase in fc2/^. For the first over- 
tones, we also found higher eigenfrequencies for the larger 
quadrupole parameter. However, the difference is much 
smaller than in the case of fundamental modes, in partic- 
ular for systems with large orbital separations. 

5. Conclusions 

We have studied the tidal effect of the companion on the 
global oscillation modes in equatorial disks around binary 
Be stars. For simplicity, we assumed the binary orbit to 
be circular and the Be disk to be inviscid, isothermal, and 
truncated at the tidal radius. Wc solved linearized equa- 
tions for global m = 1 perturbations in a three-dimensional 
Be disk with a power-law density distribution, with the 
rigid wall inner boundary condition, which is applicable 
to systems where material is ejected continuously, so there 
is no gap between the star and the disk. 

We have obtained progradc fundamental modes even 
when the quadrapole contribution to the potential is neg- 
ligible. This confirms the results of Ogilvie (2008). In 
our study of binary Be star systems, the modes are well 
confined when the disk is larger than a few tens of stellar 
radii. 

We have found that the oscillation period increases 



with increasing binary separation and/or decreasing bi- 
nary mass ratio. The effect is, however, only appreciable 
for small binary separation. In systems with large binary 
separations, the fundamental mode is well confined to the 
inner part of the disk, so the eigenfrequency no longer 
depends on the binary parameters. 

Be stars sometimes show a cavity between the star and 
the disk when the disk is being dissipated. In order to 
study the global oscillation modes in such disks, we have 
solved the perturbation equations with a free inner bound- 
ary condition, without changing the density distribution. 
We have obtained much higher eigenfrequencies with the 
free inner bounday condition than with the rigid one, as 
in Papaloizou & Savonije (2006) and Ogilvie (2008). The 
result implies an interesting possibility that the V/R vari- 
ability time-scale is shorter when the disk is dissipating 
than when it is forming or persistent. 

In this paper we have assumed a power-law density dis- 
tribution. In the disk dissipation stage, however, it is 
likely that the density distribution is far from a power- 
law form and is highly peaked near the inner radius. In 
a subsequent paper, we will study the effect of density 
distribution on the global m = 1 oscillations. 

FO thanks Masayuki Fujimoto for helpful discussions. 
She also acknowledges the scholarship from Ministry of 
Education, Culture, Sports, Science and Technology. We 
thank Gordon Ogilvie for kindly showing us his latest re- 
sult on the three-dimensional effect on the m = 1 mode 
confinement. 
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